R Notebook

fer$Treatment[fer$Treatment=="Anoxic to oxic"&fer$Sample_day<=13]<-"Anoxic"
fer$Treatment[fer$Treatment=="Oxic to anoxic"&fer$Sample_day<=13]<-"Oxic"

f = fer%>%
  filter(!is.na(Treatment))%>%
  ggplot(aes(x = Sample_day, y = NPOC, color = Treatment))+
  geom_point()+
  ylab("NPOC (mg/L)")+
  geom_vline(xintercept = 13.5)+
  scale_color_manual(values = c("#175676","#D58936","#A44200","#4BA3C3")) +
  facet_grid(Status~filt)+
  theme_bw()

fer%>%
  filter(!is.na(Treatment))%>%
  ggplot(aes(x = Sample_day, y = TNb, color = Treatment))+
  geom_point()+
  ylab("Total bound nitrogen (mg/L)")+
  xlab("Day of experiment")+
  geom_vline(xintercept = 13.5)+
  scale_color_manual(values = c("#175676","#D58936","#A44200","#4BA3C3")) +
  facet_grid(Status~filt)+
  theme_bw()

library(plotly)
ggplotly(f)

fer%>%
  filter(Sample_day == 9,
         filt == "Filtered",
         Treatment == "Oxic")

fer

f